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ABSTRACT 

Due to gravitational instability, an initially Gaussian density field develops 
non-Gaussian features as the Universe evolves. The most prominent non- 
Gaussian features are massive haloes, visible as clusters of galaxies. The dis- 
tortion of high-rcdshift galaxy images due to the tidal gravitational field of the 
large-scale matter distribution, called cosmic shear, can be used to investigate 
the statistical properties of the LSS. In particular, non-Gaussian properties of 
the LSS will lead to a non-Gaussian distribution of cosmic-shear statistics. The 
aperture mass (Afap) statistics, recently introduced as a measure for cosmic 
shear, is particularly well suited for measuring these non-Gaussian properties. 
In this paper we calculate the highly non-Gaussian tail of the aperture mass 
probability distribution, assuming Press-Schechter theory for the halo abun- 
dance and the 'universal' density profile of haloes as obtained from numerical 
simulations. We find that for values of Map much larger than its dispersion, 
this probability distribution is closely approximated by an exponential, rather 
than a Gaussian. We determine the amplitude and shape of this exponential 
for various cosmological models and aperture sizes, and show that wide-field 
imaging surveys can be used to distinguish between some of the currently most 
popular cosmogonies. Our study here is complementary to earlier cosmic-shear 
investigations which focussed more on two-point statistical properties. 

Key words: galaxies: clusters: general - cosmology: theory - dark matter - 
gravitational lensing 



1 INTRODUCTION 

The deflection of liglit due to tlie gravitational field of matter inhomogeneities is observable 
through the distortion of images of background galaxies. In the case of weak gravitational 
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fields, i.e., in the absence of strong gravitational lensing effects like giant arcs, the images of 
a population of background sources with known intrinsic ellipticity distribution can be used 
to statistically investigate weak gravitational lensing effects by measuring a net ellipticity. 

Cosmic shear - the line-of-sight integrated tidal gravitational field - reflects the statistical 
properties of the density fluctuation field (Gunn 1967, Blandford & Jaroszynski 1981). A 
quantitative description of this connection is given by the two-point correlation function of 
galaxy-image ellipticities, or by the rms-ellipticity within a (circular) aperture, which was 
investigated by several authors using the linear and nonlinear power spectrum of density 
fluctuations in different cosmologies (see the recent review by Mellier 1998, and references 
therein). 

Schneider et al. (1998; hereafter SvWJK) investigated cosmic shear using the Map- 
statistics which is a spatially filtered version of the projected density field. They computed 
the dispersion of Map for the linear and nonlinear power spectrum of density fluctuations 
in different cosmologies for a broad range of filter scales and showed that this quantity is 
a sensitive and 'local' measure of the power spectrum; in fact, as shown by Bartelmann 
& Schneider (1999), the mean-squared value of Map on an angular scale ^ is a very good 
approximation to the power spectrum of the projected density field at a wavenumber 
s ~ 4.25/^^. SvWJK calculated the skewness of Map in the frame of quasi-linear theory of 
structure growth and found that it is a sensitive indicator for the density parameter Qq, 
independent of the normalization of the power spectrum (see also Bernardeau, van Waer- 
beke & Mellier 1997; van Waerbeke, Bernardeau & Mellier 1999). The skewness measures 
the non-Gaussianity of the projected density field and indicates, according to its sign, an 
asymmetric positive or negative tail of the probability distribution function (PDF). The 
skewness of Map, which reflects the skewness of the three-dismensional density fluctuations, 
was found to be positive, indicating an extended tail of the PDF towards high positive values 
of Map. This behaviour is expected: the initial density contrast 6, which is assumed to be 
a Gaussian random field, becomes non-Gaussian during its evolution through gravitational 
collapse. That means the PDF of 6 attains a cut-off near 5 = — 1 and a broad positive tail, 
according to the occurrence of underdense and overdense regions. Since Map is a linear func- 
tion of the projected density field, the PDF of Map is closely related to that of 6. Therefore 
the PDF of Map reflects the nonlinear evolution of 6. 

The projected density field k in general is defined as a projection of the three-dimensional 
density contrast 6, weighted by a redshift-dependent factor which accounts for the lensing 
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geometry and the redshift distribution of the sources. The highest peaks of k are expected 
to arise from physical objects with high three-dimensional density, i.e., collapsed haloes. In 
Kruse & Schneider (1999; hereafter KS99) we have calculated the number density of haloes 
which yield a value of Map larger than a certain threshold. We found that for all cosmologies 
considered, the number density of haloes above a threshold corresponding to a signal-to- 
noise of 5 exceeds ten per square degree, for currently feasible deep optical imaging, so that 
a wide-field deep optical survey should detect hundreds of such peaks, which can then be 
used to define a mass-selected sample of dark matter haloes (Schneider 1996). 

In this paper we continue the investigation of statistical properties of the aperture mass 
by computing its cumulative probability distribution function (CPDF) for large positive 
values of Map, assuming that all of these are caused by collapsed structures. Describing the 
density profile of the haloes by the universal mass profile found by Navarro, Frenk & White 
(1996, 1997; hereafter NFW), we can assign an angular cross-section to each halo. If this 
cross-section is integrated over the abundance of haloes as obtained from Press-Schechter 
theory, one can determine the probability to measure a value of Map larger than some 
threshold, i.e., the CPDF for the tail of Map. Similar to the number density calculated in 
KS99, the amplitude and shape of this tail reflects the abundance of dark- matter haloes, 
which can be used powerful cosmological probe. 

The rest of this paper is organized as follows: in Sect. 2 we briefly review the concept of 
the aperture mass, as applied to the NFW- profile. The cross-sectional RrGcl clS db function 
of Map, and the CPDF, are derived in Sect. 3, and results are given in Sect. 5. The number 
density of peaks in the two-dimensional distribution of Map will be strongly affected, though 
in a controllable way, by noise, mainly coming from the intrinsic ellipticity distribution of 
galaxies. We consider an alternative observable for the abundance of peaks in Map in Sect. 
4. Finally, we summarize and discuss our results in Sect. 6. 

2 FORMALISM 

Following Schneider (1996), we define the spatially filtered mass inside a circular aperture 
of angular radius 6 around the point ^, 



where the continuous weight function [/{i}) vanishes for > 6. li U{'d) is a compensated 
filter function. 




(1) 
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f dT? d U{d) = 0, (2) 
Jo 

one can express Map in terms of the tangential shear 7t(^; C) position ^ + C relative to 
the point C 

Map(C) = /d2e7t(^;C) Qd^l), (3) 

where 

7t(^;C) = -Re(7(^ + C)e-'^^) (4) 

is the tangential component of the shear at relative position ^ = cos 0,^ sin 0). The 
function Q is related to U by 

Qi^) = 4 r - f/(^). (5) 

JO 

We use a filter function from the family given in SvWJK, specifically we choose that with 
/ = 1. Then writing U{^) = u{^/e)/e^, and Q{^) = q{^/e)/e^, 



u{x) = -{l-x')[^^-x'j, (6) 
and 

q{x) = -x\l-x^), (7) 

TT 

with u{x) = and q{x) = for x > 1. We will describe the mass density of dark matter 
haloes with the universal density profile introduced by NFW, 

MG^^ ' n{z) r/r,{l + r/r,Y' ^> 

with 

fid and fiv denote the present-day density parameters in dust and in vacuum energy, re- 
spectively. Haloes identified at redshift z with mass M are described by the characteristic 
density 6^ and the scaling radius = r2oo/c where c is the concentration parameter (which 
is a function of 6c) and r2oo is the virial radius, defined such that a sphere with radius r2oo 
has a mean interior density of 200 pcrit and contains the halo mass M2oo- We compute the 
parameters which specify the density profile according to the description in NFW using the 
fitting formulae given there. The surface mass density of the NFW-profile is given by (see 
Bartelmann 1996) 
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with 

fix) 



1 



— 1 



2 , 1-x ^ 

, 1 , arctanh \ , for x < 1 

X ^ Vl + x 
. 2 /x - 1 . 
1 , arctan W , tor x > 1 

and 6s = Vg/Dd- -Dd is the angular diameter distance to the lens. Introducing the critical 
surface mass density 
D 

= fl^, (12) 

with Ds and D^s being the angular diameter distances to the source and that from the 
lens to the source, we can define the dimensionless surface density (convergence) which is a 
function of source redshift 



_ S(?9) 

with 



<^,Zd,Zs) = ^ = Kofij\, (13) 



no = Hl + z) (14) 
We assume a normalized source redshift distribution of the form 

P.iz) = — ^ eM-[z/zof), (15) 



(see Brainerd et al. 1996). The mean redshift of this distribution is proportional to zq and 
depends on the parameter f3 which describes how quickly the distribution falls off towards 
higher redshifts. We will use the values /3 = 1.5 and zo = 1- For these values the mean redshift 
{z) is given by {z) = 1.505 Zq. With the distribution (|15|) we define a source distance-averaged 
surface density 

K,{i9,Zd) = f dzs Pzizs) K{i?,Zd,Zs) = Rof ( ^ ) , (16) 



with kq = J dzs Pz{zs) kq. Inserting ([l6|) in we get, after introducing polar coordinates 
for i9 and setting ^ = (C?0), which is the radial distance of the aperture centre from the 
halo center, 

Map(C) = / d^^ 



d(f) U{J^^ + C~ 2^9Ccos(20) , (17) 
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Figure 1. The aperture mass, as defined in (p^), as a function of the radial distance C, computed for five different cosmologies 
as indicated by the line types. The numbers in parentheses are the normalization erg annd the shape parameter V. The halo 
mass IS W'^^Mq/H and its redshift = 0.3. The filter radius is = 2 arcmin. 

where 0m = min (vr, arccos ^Sc~^ ))' "^^^ filter function chosen above , the 0-integration 
can be performed analytically. 

In Figure |l| we plot the aperture mass of a halo with mass M = IO^^Mq/H at redshift 
Zd = 0.3 as a function of the radial distance (, using a filter scale of ^ = 2', for five 
cosmological models. If we fix the parameters (M, z^, 0), the aperture mass is a monotonically 
decreasing function of ( between ( = and the first root of Map. 



3 THE TAIL OF Map-STATISTICS 

In this section we calculate the CPDF P(> Map, 9), i.e., the probability to find an aperture 
mass larger than Map using a filter with radius 6. We concentrate on values of Map which are 
much larger than the rms value of Map, i.e., we consider only the far tail of the probability 
distribution. We may assume that such high values of Map are caused exclusively by collapsed 
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haloes. Thus, from the properties of such haloes, together with their abundance, we can then 
determine the CPDF. 

Assuming a halo characterized by its mass M and redshift Zd, we can invert the function 
(p!7|) for a given value of Map and a fixed filter radius 9 (see Fig. |T]). As a result we get the 
separation ( = ({Mi^p,9, M, z^). Owing to the monotonic behaviour of Map between C = 
and its first root, separations smaller than the one obtained by inversion correspond to 
aperture masses larger than the threshold Map. Therefore, ( defines for each halo an angular 
cross section 

c7(Map, 9, M, Zd) := vr ('(Map, 9, M, Zd), (18) 

which represents the halo target area for detecting a weak lensing signal with an aperture 
mass larger than Map. The cross section (|18]) is non-zero if the aperture mass measured in 
the halo centre is larger than the threshold Map (see Fig. |l|). 

The CPDF is now obtained by summing up the cross sections of all haloes within a unit 
solid angle; this yields 

Pi> M^,9) = ^ fdzd ^^-^ Dlizd) 

xJdM A^haio(M, Zd) a(Map, 9, M, Zd), (19) 
with 



E{zd) = Vfid(l + Zd)^ + (1 - fid - fiv)(l + Zd)^ + Ov. (20) 

Nhaio{M,Zd) dM is the comoving number density of haloes with mass within dM about 
M at redshift Zd- We assume that A^haio(^5 -^d) can be obtained from Press and Schechter 
(1974) theory (see, e.g. Lacey & Cole 1993, 1994). The M-integral in (|T^) extends from a 
lower threshold Mi to infinity, where a vanishes for M < Mt (see KS99). 



4 NUMBER DENSITY OF HALOES, SELECTED BY Map AND SIZE 

In KS99, we have calculated the number density of haloes A^(> Map, 9) with an aperture 
mass > Map, using the same physical model as described above. Considering Map(^^) as 
a two-dimensional map, A^(> Map, 9) yields, at fixed filter scale 9, the number density of 
peaks in this map with amplitude > Map. Such peaks may be generated by noise where the 
major contribution comes from the intrinsic ellipticity distribution of galaxies. Intuitively, 
one might expect that high peaks are less effected by noise than smaller ones; this motivates 
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Figure 2. The cross section radius obtained from inverting jl^ ) as a function of the lens mass for the same cosmological models 
as in Figure |^. ( p^ ) is inverted for the aperture mass Map = 0.04, a filter radius 9 = 2 arcmin and a halo redshift = 0.3. 
The cross section of haloes characterized by the fixed parameters written in the panel is zero for halo masses smaller than 
Mt ~ 10^3.8^^^^/;^ 

US to consider the number density of haloes with aperture mass > Map and cross-sectional 
area > a = 7i(^, which is again obtained from summing over all haloes per unit solid angle, 

N{> Map, > Ct, ^) = J- / d.d Dliz,) 

X JdM iVhalo(M, ^d) H[C(Map, e, M, Zd) - Ct], (21) 

where H is the Heaviside step function. The integrand is non-zero only for M > Mt, where 
Mt = Mt(Ct, 2;d, Map, 6') is the mass obtained by inversion of the function shown in Fig. |^. 
For Ct = 0, N{> Map, > Ct, 0) = N{> Map, 9). The value of Ct for which N{> Map, > Ct, 0) 
decreases to about 1/2 of A^(> Map, 6') yields the characteristic size of peaks corresponding 
to a given Map; this value is expected to be < 6. E.g., for the EdS(0.6,0.25) cosmology we 
obtain A^(> 0.04, > 0.45', 2') = 4.6 which can be compared to A^(> 0.04, 2') = 9.4 computed 
in KS99. For this example the characteristic size of a peak is roughly 1/4 of the filter radius. 
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5 RESULTS 

In this section we consider P{> Map, 6) for different cosmological models in the regime 
where the PDF of Map describes the non-linear evolution of the density field. Furthermore, 
we use the observable (^TJ) to get additional constraints for the cosmological parameters. We 
perform our calculations for the same five cosmological models as shown in Figure |^. For 
three of them, the power spectrum is approximately cluster normalized, which corresponds to 
as ~ 0.6 for an Einstein-de Sitter universe (EdS, Qd = 1, = 0) and as = 1 for both an open 
universe (OCDM, = 0.3, = 0) and a spatially flat universe with cosmological constant 
(ACDM, Qd = 0.3, fiv = 0.7). For these models we use the shape parameter F = 0.25 which 
yields the best fit to the observed two-point correlation function of galaxies (Efstathiou 
1996). The remaining two EdS models have higher normalization {as = 1, approximately 
corresponding to COBE normalization) or a different shape parameter (F = 0.5). 
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In the absence of lensing, the expectation value of Map vanishes, and its dispersion 
depends on the intrinsic elhpticity distibution, as calculated in Schneider (1996), 



P(> Map, e) = po exp 



(24) 



We take n = 30 arcmin"^ and = 0.2 for the number density of background galaxies and 
their intrinsic ellipticity distribution, respectively, as representative values in the following. 
The signal-to-noise ratio of Map is then defined as 

SM ^ (23) 

In Fig. 1^ we plot the probability (^) for four different filter scales for the cosmologies 
introduced in the beginning of this section, for values of Map between Mq and 2Mq, where 
the value of Mq was chosen to correspond to a signal-to- noise ratio of five, Mq = 5ac{6). 
According to this figure the probability P{> Map, 6) can be well approximated by an expo- 
nential, 

(Map - Mq) 
C 

where po and c are fit parameters. From (^4]) it is clear that po is the probability to find a value 
of Map larger than Mq. We determined Pq and c simply by assuming that the logarithm of 
P(> Map, 6) in the interval [Mq, 2Mq\ follows a straight line. The maximal relative deviation 
between the probabihty (0) and its approximation (|^) in Fig. |^ is about 3 percent. 

If we compute the probability according to the approximation (p^, all information on 
cosmology is contained in the parameters po and c. If we attempt to constrain the cosmolog- 
ical parameter set, we have to compute dispersions for the fit parameters. In the following 
we employ a maximum likelihood analysis to derive confidence levels in the parameter space 
{c,Po}- 

In order to obtain a likelihood function we have to specify a probability distribution for 
finding a particular set of values {Map}, i G [1, A^f] of N{ statistically independent aperture 
mass values. As an example, we shall assume that we have an image with side length L. 
SvWJK showed that about A'f = {L/{26)Y statistically independent fields can be placed on 
the image. This means that we consider two fields to be statistically independent if their 
angular separation is about two filter radii since then the correlation coefficient of the two 
fields is below 1 percent (see SvWJK). We suppose that this sample contains A^> values of 
Map which are above the threshold Mq. This subsample consists of aperture mass values 
M^^p > Mq, j G [1,A^>]. Now we may ask for the probability to find A^> G [l,A^f] values 

© 1998 RAS, MNRAS 000, §-0 



The non-Gaussian tail of cosmic-shear statistics 11 



0.0014 

0.0012 

0.0010 

^ 0.0008 

0.0006 

0.0004 

0.0002 
0.0000 



1 1 1 1 1 1 1 1 1 


1 1 1 1 1 1 1 1 1 1 1 1 1 
OCDM (1.0,0.25) 


L 0-1' 


( ^ 




: ACDM (1.0,0.25) 






/• , 

I ^, 
s, 

\. 




') ^ 
y : 


: EdS (0.6,0.5) ^ 








+ z:> 








EdS (0.6,0.25) ■ 



0.008 
0.006 - 

0.004 - 

0.002 - 
0.000 1. 



© = 2' 



0.005 



0.015 



0.025 



0.006 0.010 0.014 O.Oli 



0.01 4 
0.01 2 : 
0.01 : 

0.008 : 
0.006 : 

0.004 - 



0.035 r 



8 = 3' 




0.006 0.008 0.010 0.012 0.014 
c 



0.030 

0.025 7 
0.020 7 

0.01 5 7 

0.01 r 
0.005 I 



1 8 = 5' 



^[^ (^< T \ 




0.0030 



0.0050 



0.0070 
c 



0.0090 



Figure 4. The fit parameters of the exponential (^) denoted by crosses computed as described in the text together with their 
l-tr error eUipses which are defined by the dispersions ( psj ) and (po|). For all panels the signal-to- noise ratio threshold Mq/ctc 
is 5. Haloes with redshift in G [0, 1] are considered. We use the source redshift distribution ( p^ ) with /3 = 1.5 and zq = 1. 
The thin and the thick curves describe a 25 deg'^ and a 100 deg^ survey, respectively. The cosmological models are indicated 
in the upper left panel. The different panels are obtained by varying the filter scale, 9 = 1, 2, 3, 5 arcmin. We do not plot the 
EdS(1.0,0.25) models because it is well separated from the other cosmologies (po is about a factor of 5 larger compared to the 
probabilities of the remaining model). 

above the threshold if we have a sample of A^f fields. For each of the A^f aperture mass 
measurements the event Map > Mq is supposed to occur with probability po. Obviously this 
random process can be described by a binomial distribution, 



P(iV>|iVf,po) 



iV> 

Po 



:i-po)^'-^>. (25) 

Of course, we can hope to estimate the two parameters po and c only for observations for 
which A^> 3> 1; in particular, if A^> = 0, the parameter c will be completely undetermined. 
Thus, we want to consider only realizations of p5|) with A^> > 1, and therefore renormalize 

P{N^\Ni,po), 



f J2 P(iv>|iVf,Po) = i, 



(26) 
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Figure 5. The top left plot is the same as in Fig. ^ for L = 10 deg. In the top right plot all parameters are unchanged but only 
halo redshifts g [0.6, 1] are considered. The bottom left plot is the same as the top right one but halo redshifts are from the 
interval £ [0,0.4]. In the bottom right plot all sources are located at redshift z^ = 3; the other parameters are the same as 
in the top left plot. As in Fig. ^ we do not plot the high normalized EdS model which again is well separated from the other 
models. In the top right panel the value of po for the EdS(0.6,0.25) model is so small that the normalization constant / differs 
from unity appreciably; we have therefore not included this model in the figure. The cosmological models are indicated in the 
upper left panel 



where 
/ 



(27) 



l-(l-po)^f' 

In addition, we consider the probabihty distribution of Map itself for values above the thresh- 
old Mq. The normalized PDF of Map > Mq can be obtained from (^) by 



p(Map) = - 



dM 



-P(> Map, I 



ap 



1 

- exp 

c 



(Map - Mo) 



Combining both probabilities, we obtain the likelihood function 
L(Mi,...,MiJ>|po,c) = P(iV>|iVf,po) n P(^ip)' 



(28) 
(29) 
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which describes the probabihty to find A^> aperture masses > Mq where those above the 
threshold follow the distribution (^). From (p9D we can derive maximum likelihood estimates 
for both fit parameters. Denoting these estimates by (c, po) we obtain by differentiating ( pQ]) 
w.r.t. c and po 

c = TrE^ip-^o (30) 

> j=l 

and 

PO = f ; (31) 

these results are of course not unexpected. In order to show that ( pOD and (0) are unbiased 
estimators we have to evaluate their ensemble average which can be performed by applying 
the operator 



7V>=1 \^>. 



X 



n /, dMipp(Mip) (X). (32) 



Mo 



Averaging (|30D and (0) with (p2D yields 

(c) = P(c) = c and (po) = P(po) = Pof- (33) 

Because of the second of eqs.(^), po is an unbiased estimator for po only if / ^ 1. In the 
following we will consider only values of Pq and Nf which guarantee that / ^ 1 is satisfied, 
otherwise the statistics is not good enough to determine the parameters accurately anyway. 
The correlation between the two estimators can be calculated from 

{cpo) = F{cpo) = cpof. (34) 

Since we have / ~ 1 these two estimators are not correlated. In the following we set / = 1. 
Dispersions of the estimators are defined by 



^Po = ^{PD - {Po? (35) 
and 

a, = v'(a^) - {d?. (36) 

The calculation of the dispersions involves the application of the operator (|3^ to the squared 
estimators. After some algebra we obtain 
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and 




(38) 



(Tc = C 



^ P(iV>|iVf,po) iV: 



(39) 



and 




iVf 



Po 



(40) 



As expected, the rms values of the estimators decrease with increasing image size. This 



compute the binomial distribution for a given Ny. The offset for the recurrence relation, 
(7(1) = A^fPo(l — Po)^'~\ rapidly decreases for large A^f and therefore the sum in (^). 

Note that the rms values of the estimators depend only on their means and iVf which is de- 
termined by the image size for a given filter scale. Given that the mean values (c, po) depend 
on the filter radius 9 and the threshold Mq we can vary the parameter triple (6', Mq, Af) 
to obtain a significant difference between the various cosmological parameters. Further- 
more, redshift information coming from possible redshift measurements of the haloes and 
the sources can be used to maximize the difference between various cosmologies. 

In Fig. ^ we plot the mean values of the fit parameters (denoted by crosses) for the filter 
radii 9 = 1, 2, 3, 5 arcmin for four cosmological models. The threshold Mq is defined by a 
signal-to-noise ratio of 5. For each cosmology and filter scale we use a 25 deg^ (thin lines) 
and 100 deg^ (thick lines) survey to compute the dispersions of the estimators (|3T|) and (|30D 
which define the l-cr error ellipses in Fig. ^. As expected, if we double the side length of 
the image the dispersions of c and po become smaller by about a factor of two. The COBE 
normalized EdS model is well separated for all filter scales. Since the probability po in this 
model is about a factor of 5 larger than that of the remaining cosmologies we do not plot this 
model. The possibility of distinguishing cosmologies strongly depends on the filter scale. For 
an aperture with 1 arcmin filter radius we can clearly distinguish the EdS(0.6,0.25) and the 
low density model without cosmological constant from the other cosmologies if we use the 
survey with the larger area. If we enlarge the filter radius, the differences between the two 
low-density models and the EdS model with large shape parameter become smaller whereas 
the EdS(0.6,0.25) model remains distinguishable from all the other cosmologies. 



behaviour is obvious for (|40|) and can be seen for (|39|) if we use a recurrence relation to 
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In Fig. ^ we demonstrate the effects on the dispersions if we change the halo redshift 
integration range and the source redshift ^s, where we assume all sources to be located at 
the same redshift. For reference, we choose the upper left plot from Fig. ^ with L = 10 
deg. If we shrink the halo redshift interval to = [0.6, 1] we obtain the upper right plot. 
Because of the stronger evolution of the halo number density for high redshifts in the other 
cosmologies, the EdS(0.6,0.5) model can now be distinguished from the low-density models. 
The EdS(0.6,0.25) model has a value po which is so small as to lead to a significant deviation 
of / - as defined in (^) - from unity, so that the expressions (|39|) and (^O]) are no longer 
valid, and therefore we have not plotted this model. 

In the bottom left plot, compared to the upper right one, we only consider haloes having 
redshifts in the interval = [0,0.4]. In this redshift range all cosmologies considered are 
distinguishable. The reason for the low-density models now being separated is the increas- 
ing difference between the rich cluster mass functions of both models if we choose smaller 
redshifts. 

In the remaining plot, in comparison with the top left one, all sources are assumed to 
have the same redshift Zs = 3 which is about twice the mean source redshift used in the 
other panels. Because of the improved efficiency of the weak lensing signal, the probability 
is increased by about a factor of 4. For this very deep survey all cosmologies are clearly 
separated. From Figs. § and ^ it becomes clear that we need large deep surveys in order to 
be able to distinguish clearly between the cosmological models considered, using only the 
statistics applied here. 

If we want to compute P{> A^ap? ^) as expected from real observations, we have to 
consider possible sources of error which may change the theoretical value of Map, because 
the cross sections are given via aperture mass measurements (see Fig. |l|). To account for 
the intrinsic ellipticity distribution of galaxies, we can proceed in the same way as for the 
observable A^(> Map,6) in KS99. Essentially this means that we convolve P{> Map,^) 
with a Gaussian defined by the dispersion of the intrinsic galaxy ellipticity distribution 
of the sources. As explained in SvWJK, this dispersion is expected to dominate the noise 
in a measurement of Map. The probability obtained from the convolution can be directly 
compared with observations. We have checked that, as expected from the results of KS99, 
the convolution only slightly enhances the theoretical values and does not change the shape 
of P(> Map, 6) appreciably in the range of Map considered here. Therefore the fit formula 
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Table 1. The number of haloes per square degree with aperture mass larger than Map = 0.04 and cross-section radius greater 
than <^t = 0.8 arcmin, as defined in (^), computed for five cosmological models. The filter radius is d = 2 arcmin. The aperture 
mass and the dispersion determined by the filter radius correspond to a signal-to-noise ratio of 5. 
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o» 


Af(> 0.8', > 0.04,2') 


1.0 





0.25 


0.6 


0.47 


1.0 





0.25 


1.0 


11.9 


1.0 





0.5 


0.6 


1.16 


0.3 





0.25 


1.0 


2.22 


0.3 


0.7 


0.25 


1.0 


1.67 



we obtained remains valid, and the maximum likelihood method can be applied to derive 
dispersions for the fit parameters. 

In Table [I| we show the number density of haloes with aperture mass larger than M^p = 
0.04 and cross-section radii exceeding the threshold Ct = 0.8 arcmin for the filter radius 
6 = 2 arcmin. According to the fact that the rich cluster mass function shows most clearly 
the cosmology dependence of the non-linear evolution of the density field and since cross- 
section radii above 0.8 arcmin correspond to the most massive non-linear objects (see Fig. 
H), A^(> Map, > 0.8', 6') better distinguishes cosmologies than N{> Ma.p,9) for the same 
Map and 6 (see KS99). A drawback when using the observable (|2T|) is the larger image size 
required for the detection of significant differences between various cosmologies. From the 
numbers in Table |l| we infer a survey area of 25 deg^ which is needed to distinguish the 
cosmologies considered here significantly, i.e., with no overlapping Poissonian error bars. 

6 DISCUSSION AND CONCLUSIONS 

In this paper we computed the highly non-Gaussian tail of the probability distribution of 
the aperture mass Map resulting from lensing by the large-scale structure. The CPDF is 
obtained by summing up the cross sections of dark matter haloes with assumed spherical 
mass density following a universal NFW density profile. We used Press-Schechter theory to 
compute the number density of massive haloes which cause the extended tail in the PDF of 

Map. 

The number density of haloes with large cross-sectional radii, or in other words, with large 
masses, is a sensitive measure for constraining the cosmological parameter set (see KS99, and 
references therein). We showed that the number density in a mass-selected sample of haloes 
with cross sectional radii above 0.8 arcmin, which corresponds to halo masses ~ 10^^Mq//i, 
is measureable in all cosmologies considered. Especially, from a deep, high-quality imaging 
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survey of 25 deg^, some of the currently most popular cosmologies can be distinguished by 
the varying number of haloes with the selected cross section threshold. The expected number 
of massive haloes varies from ~ 12 for the EdS(0.6,0.25) to ~ 300 for the high normalised 
EdS model if we use the above-mentioned survey. 

The main result of this work is an analytical formula describing the PDF of the aperture 
mass, which turns out to be closely approximated by an exponential for values of Map well 
above its rms. This fact allows a simple maximum-likelihood analysis for distinguishing 
various cosmological models, using the non-Gaussian tail of Map only, as demonstrated in 
Sect. 5. 

Modifications of the present investigation can well be imagined and may turn out to yield 
fruitful results in future. One is the use of (photometric) redshift estimates which allow a 
more precise measurement of the shear, owing to the redshift-dependent geometrical factors 
entering the projected surface mass density. If the haloes found by our method turn out 
to be associated with galaxy concentrations, their redshift can be estimated, and the halo 
abundance as a function of Map and redshift can be obtained, allowing to greatly refine the 
statistics considered here and in KS99. 

Our approach to attempt measuring cosmological parameters is complementary to inves- 
tigations employing two-point statistical measures, such as the rms shear in (circular) aper- 
tures or the shear two-point correlation function (e.g., Blandford et al. 1991; Kaiser 1992, 
1998; Villumsen 1996; Jain & Seljak 1997; Bartelmann & Schneider 1999), and those which 
consider the skewness of cosmic shear statistics as a powerful handle on (e.g., Bernardeau 
et al. 1997; van Waerbeke et al. 1999; Jain et al. 1999). In contrast to these lower-order sta- 
tistical measures, the highly non-Gaussian features of cosmic shear are expected to be less 
affected by noise, whose main contribution is the intrinsic ellipticity dispersion of source 
galaxies. 

We note that the cosmology dependence of our results are mainly (but not exclusively) 
through the abundance of dark matter haloes as a function of mass and redshift. Thus, our 
statistics provides a fairly direct measure of this cluster abundance. In a future work, we 
shall attempt to relate Map to the true three-dimensional mass of haloes as determined in 
numerical A^-body simulations, and thus relate the aperture mass statistics directly to the 
mass function of haloes. 

It should be stressed that, whereas the aperture mass is most likely not the optimal 
statistics to measure the cosmic shear power spectrum (see Kaiser 1998, and Seljak 1998, for 
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different approaches), it is a particularly convenient measure for highly non-Gaussian, spa- 
tially localized features which can be obtained locally and directly from the observed image 
ellipticities and whose noise properties can be straightforwardly investigated. In particular, 
cosmic shear measurement using Map, and the search for mass-selected haloes (Schneider 
1996; KS99) are just two aspects of the same underlying physics and statistics. 

The validity of the various approximations which enter KS99 and the current study needs 
to be investigated in more detail. In a forthcoming paper (Reblinsky et al., in preparation) 
we will apply the aperture mass statistics to the same numerical simulations used in Jain 
et al. (1999), which combine very large N-body simulations with ray-tracing methods. Such 
numerical simulations are indispensible not only to assess the accuracy of the analytical 
approximation, but also to study statistical estimators which cannot be calculated analyti- 
cally. Given the highly non-Gaussian nature of the projected density field resulting from the 
evolved LSS density distribution, it is by no means clear how to optimally and robustly dis- 
tinguish between different cosmological models. The use of the far tail of the Map-statistics 
should be viewed as one of several useful tools. 
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